Kramers degeneracy in a magnetic field and Zeeman spin-orbit coupling 

in antiferromagnetic conductors 
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In this article, I study magnetic response of electron wavefunctions in a commensurate collinear 
antiferromagnet. I show that, at a special set of momenta, hidden anti-unitary symmetry protects 
Kramers degeneracy of Bloch eigenstates against a magnetic field, pointing transversely to staggered 
magnetization. Hence a substantial momentum dependence of the transverse g-factor in the Zeeman 
term, turning the latter into a spin-orbit coupling, that may be present in materials from chromium 
to borocarbides, cuprates, pnictides, as well as organic and heavy fermion conductors. 

PACS numbers: 75.50.Ee 



Antiferromagnetism is a frequent occurrence in mate- 
rials with interesting electron properties: it is found in 
elemental [l[ and binary Q solids, in numerous borocar- 
bides 0] , in doped insulators such as cuprates Q , in iron 
pnictides [j| , in various organic [|| and heavy fermion Q 
compounds. The physics of an antiferromagnetic state in 
these materials has been a subject of much research. 

In this article, I study response of electron states in an 
antiferromagnet to a weak magnetic field. I concentrate 
on the simplest case: a commensurate collinear antiferro- 
magnet, shown schematically in Fig. [1] where magneti- 
sation density at any point is parallel or antiparallel to a 
single fixed direction n of staggered magnetisation, and 
changes sign upon primitive translation of the underlying 
lattice. 

In a paramagnet, double degeneracy of single-electron 
eigenstates is usually attributed to symmetry under time 
reversal 9 - and, indeed, perturbations that break (such 
as ferromagnetism or magnetic field) tend to remove the 
degeneracy. Yet mere violation of does not preclude de- 
generacy: in a commensurate centrosymmetric Neel an- 
tiferromagnet, as in a paramagnet, all Bloch eigenstates 
enjoy Kramers degeneracy Q in spite of time reversal 
being broken in the former, but not in the latter. 

In an antiferromagnet, staggered magnetisation sets 
a special direction n in electron spin space, making it 
anisotropic. Magnetic field along n removes the degen- 
eracy, as it does in a paramagnet. By contrast, in a 
transverse field, the symmetry remains high enough to 
protect Kramers degeneracy at a special set of momenta. 
Generally, in d spatial dimensions, full degeneracy man- 
ifold is (d — l)-dimensional; at its subset, degeneracy is 
dictated by symmetry. This is in marked contrast with 
a paramagnet, where arbitrary magnetic field lifts the 
degeneracy of all Bloch eigenstates. 

I show that, at a subset of the degeneracy manifold 
above, it is a hidden anti-unitary symmetry that protects 
Kramers degeneracy of Bloch states in a transverse field. 
The degeneracy gives rise to a peculiar spin-orbit cou- 
pling, whose emergence and basic properties, along with 
the degeneracy itself, are the main result of this work. 



Kramers degeneracy of special Bloch states in a trans- 
verse field means, that the transverse component g± of 
the electron g-tensor vanishes for such states. Not being 
identically equal to zero, g± must, therefore, carry sub- 
stantial momentum dependence, and the Zeeman cou- 
pling Tizso must take the form 



Uzso = ~HB [g\ (H|| • a) + g± (p) (H ± • a)] 



(1) 



where H|| = (H • n)n and = H H | arc the longi- 
tudinal and the transverse components of the magnetic 
field with respect to unit vector n of staggered magneti- 
sation, fis is the Bohr magneton, while g\\ and g±(p) are 
the longitudinal and the transverse components of the 
g-tensor. 

This very momentum dependence of <?j_(p) turns the 
common Zeeman coupling into a Zeeman spin-orbit inter- 
action Jizso ©j whose appearance and key properties 
are at the focus of this work. Zeeman spin-orbit coupling 
may manifest itself spectacularly in a number of ways, 
which will be mentioned below and discussed in detail 
elsewhere. 

Symmetry properties of wave functions in magnetic 
crystals have been studied by Dimmock and Wheeler 0] , 
who pointed out, among other things, that magnetism 
not only lifts degeneracies by obviously lowering the sym- 
metry, but also may introduce new ones. This may hap- 
pen at the magnetic Brillouin zone (MBZ) boundary, un- 
der the necessary condition that the magnetic unit cell 
be larger, than the paramagnetic one 0]. 

For a Neel antiferromagnet on a square lattice, re- 
sponse of electron states to magnetic field was studied 
in [10j by symmetry analysis, and in ll| within a weak 
coupling model. The present work revisits (Io| . extends 
it to an arbitrary crystal symmetry and shows, that the 
picture is more rich than envisaged by the authors. At 
the same time, the present work extends Q by allow- 
ing for external magnetic field - to show how, at spe- 
cial momenta, Kramers degeneracy may persist even in 
a transverse magnetic field. 

Antiferromagnetic order couples to the electron spin er 
via exchange term (A r • er), where A r is proportional to 
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the average magnetisation density at point r. Nonzero 
A r changes sign under time reversal 9, and removes the 
symmetry under primitive translations T a , thus reducing 
the symmetry with respect to that of paramagnetic state. 
In a doubly commensurate collincar antiferromagnet, A r 
changes sign upon T a : A r+a = — A r , while T a leaves 
A r intact: A r+ 2 a = A r . Even though neither 9 nor T a 
remain a symmetry their product 9T a does (see Fig. Q]). 
In a system with inversion center, so does 9T a X, where 
X is inversion. 




FIG. 1: (color online). Doubly commensurate collinear anti- 
ferromagnet on a simple rectangular lattice. In the absence 
of magnetism, time reversal and primitive translations 
and T y , shown by dashed arrows, are symmetry operations. 
In the antiferromagnetic state, neither of the three remains 
a symmetry, but the products 8T X and 0T V , shown by solid 
arrows, do, as illustrated by filled spin arrows. Small dashed 
rectangle at the center is the Wigner-Seitz cell boundary in 
the paramagnetic state, while the shaded hexagon is its an- 
tiferromagnetic counterpart. Notice that neither of the point 
group operations interchanges the two sublattices, hence any 
point symmetry of the lattice, including inversion T, remains 
a symmetry of the antiferromagnetic state. 



Combined anti-unitary symmetry 9T a X induces 
Kramers degeneracy [8]: If |p) is a Bloch eigenstate 
at momentum p, then 9T a X\p) is degenerate with |p). 
Since 9 and X both invert the momentum, both |p) and 
6>T a Z|p) carry the same momentum label p. At the 
same time, |p) and 6>T a Z|p) are orthogonal: recalling 
that (T a I) 2 = -9' 2 = 1, one finds Q 



(p\9T a X\p) - -(p\9T a X\p). 



(2) 



Thus, in spite of broken time reversal symmetry, in a 
centrosymmetric commensurate Neel antiferromagnet all 
Bloch states retain Kramers degeneracy. 

Generally, magnetic field H lifts this degeneracy. How- 
ever, in a purely transverse field, hidden anti-unitary 
symmetry may protect the degeneracy at a special set 
of points in the Brillouin zone, as I show below. 



In a commensurate collinear antiferromagnet in mag- 
netic field H, electron Hamiltonian has the form 



H = Ho + (A r ■ a) - (H • <x), 



(3) 



where 'paramagnetic' part TCq is invariant under indepen- 
dent action of T a and 9, and g\iB is set to unity. In the 
absence of the field, all Bloch eigenstates of Hamiltonian 
(|3"|) enjoy Kramers degeneracy by virtue of Eqn. @. 
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FIG. 2: (color online). Relative orientation of A r , A r + a , Hy 
and Hi. Notice that 8 flips both A r and H, while T a leaves 
H intact, but inverts A r . 

Consider symmetries of Hamiltonian ((3|), involving a 
combination of an elementary translation T a , time re- 
versal 9, or a spin rotation U m (<^) around axis m by 
angle <j>. The relative orientation of A r , Hy and is 
shown in Fig. [2] 

Transverse field H j_ breaks the symmetries U n (</>) and 
T a 9 (both change Hi), but preserves U„(7r)#T a , their 
combination at 4> = n. Acting on the exact Bloch state 
|p) at momentum p, this combined anti- unitary opera- 
tor creates a degenerate partner eigenstate U n (7r)6'T a |p), 
which is orthogonal to |p) everywhere in the Brillouin 
zone, unless p belongs to a paramagnetic Brillouin zone 
boundary: 



(p|U n (7r)0T a |p> = e- 2l P- a (p|U n (7r)0T a |p 



(4) 



Equation ([4]) can be derived analogously to Eqn. ([2]) as 
soon as one observes, that [U n (7r)0T a ]" = T a 2 . 

Notice, however, that eigenstate U n (7r)#T a |p) carries 
momentum label — p rather than p. Generally, these two 
momenta are different. However, an important exception 
takes place at the magnetic Brillouin zone boundary, if 
there is a unitary symmetry U, transforming — p into 
a momentum, equivalent to p modulo reciprocal lattice 
vector Q of the antiferromagnetic state Q : 



Up = p + Q. 



(5) 



In this case, eigenstate WU n (7r)0T a |p) carries momen- 
tum label p + Q = p, is degenerate with |p) and orthog- 
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onal to it, thus explicitly demonstrating Kramers degen- 
eracy at momentum p in a transverse field. In the sim- 
plest case, as in Fig. [3] below, U is the unity operator. 

Additional insight into the locus of states, whose de- 
generacy persists in a transverse magnetic field, is af- 
forded by weak-coupling Hamiltonian of a single electron 
in a doubly commensurate collinear antiferromagnet. Let 
Q be the antiferromagnetic ordering wave vector (see the 
examples below); A r creates a matrix element (A • a) 
between the Bloch states at momenta p and p + Q; for 
simplicity, I neglect its possible dependence on p. In 
magnetic field H, and at weak coupling, Hamiltonian (J3]) 



takes the form 



H 



11] 



e p - (H ■ a) 



(A- a) 



(A • a) e P+Q — (H • er) 



(6) 



where e p and c p +q are single-particle energies of TLq in 
(J3j) at momenta p and p + Q. 

In a purely transverse field H^, the spectrum of this 
Hamiltonian is simply 



£ p = t 7p ± v /A 2 + [C p T(H ± - ( t) 



where r) p 



he p+Q 



and C p 



(7) 

. ,v.... M :> - . ■ Equation ([7]) 
illustrates several points. Firstly, at half-filling, a gap 
of size 2A opens at the chemical potential. Secondly, in 
the absence of magnetic field, each eigenstate is indeed 
doubly degenerate, in agreement with the arguments, en- 
capsulated in Eqn. ([2]). Finally, Eqn. (JT]) shows, that the 
degeneracy persists in a transverse field (and, therefore, 
<7_i_(p) in Eqn. |l| vanishes) whenever £ p = 0. Barring a 
special situation, this equation defines a surface in three 
dimensions, a line in two, and a set of points in one. Fur- 
thermore, as shown above, this manifold must contain all 
points, satisfying Eqn. ([5]). 

Notice that transverse magnetic field not only intro- 
duces the last term in Hamiltonian ©, but also tilts the 
sublattices. However, the resulting magnetisation has the 
same symmetry as the field, and thus does not remove 
the degeneracy. 

Consider examples. In one dimension, magnetic Bril- 
louin zone boundary reduces to two points p = 
which in fact coincide modulo antiferromagnetic wave 
vector Q = £ that is also the reciprocal lattice vector of 
the antiferromagnetic state (see Fig. [3]). In terms of the 
general condition ([5|), this is the simplest case: U = 1. 

As a result, at p = ±f-, the two exact Bloch states in 
a transverse field, |p) and #T a U n (7r)|p), correspond to 
the same momentum p, and are degenerate by virtue of 
#T a U n (7r) being a symmetry. Equation ([4]) guarantees 
their orthogonality, thus protecting Kramers degeneracy 
at momentum p = against transverse magnetic field. 

Now, consider a two-dimensional antiferromagnet of 
simple rectangular symmetry, with the ordering wave 
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FIG. 3: (color online). The paramagnetic (p = ±-), and 
the antiferromagnetic (p = Brillouin zone boundaries 

of a one-dimensional Neel antiferromagnet. In the antiferro- 
magnetic state, the two points p = are identical modulo 
the antiferromagnetic reciprocal lattice vector Q = — . At 
these two points, anti-unitary symmetry U n (7r)T a # protects 
Kramers degeneracy against transverse magnetic field. 



vector Q = (tt,tt), as shown in FigfT] In a transverse 
magnetic field, degeneracy persists at a line in the Bril- 
louin zone, by virtue of Eqn. (|7|). I will show that, in 
the rectangular case, the degeneracy line must contain 
point E (i.e. the star of point p = Q/2) at the center 
of the magnetic Brillouin zone boundary (see Fig. [Ua)). 
Consider a Bloch state |p) at momentum p in a trans- 
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FIG. 4: (color online). Geometry of the problem, (a) The 
Brillouin zone for a simple rectangular lattice (the rectan- 
gle), and its antiferromagnetic counterpart (MBZ, shaded 
hexagon). Thick curve (red online), passing through point 
E, shows a typical degeneracy line g_i_(p) = 0. At the MBZ 
boundary, only momentum p at point E is equivalent to p 
modulo the reciprocal lattice vector of the antiferromagnetic 
state. For a generic p', shown by the dashed arrow, this is 
no longer true, (b) The Brillouin zone of a simple square lat- 
tice and its antiferromagnetic counterpart (shaded diagonal 
square). The degeneracy line must contain the entire MBZ 
boundary, shown in red online. 

verse field. As discussed above, eigenstate #T a U n (7r)|p) 
at momentum — p is degenerate with |p) and, according 
to Eqn. (HJ), must be orthogonal to it unless (p • a) is an 
integer multiple of ir. At points E (the star of p = Q/2), 
X, and Y, momenta p and — p coincide modulo a recip- 
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rocal lattice vector of the antiferromagnetic state. How- 
ever, at points X and Y (as well as at the entire vertical 
segment of the MBZ boundary in Fig. Ufa)), (p • a) is an 
integer multiple of 7r; hence |p) and 6*T a U n (7r)|p) are not 
obliged to be orthogonal there as per Eqn. (g]) . Thus, E 
is the only point at the MBZ boundary where the two 
degenerate states |p) and #T a U n (7r)|p) are orthogonal 
and correspond to the same momentum. Dashed arrows 
in figure |4ja) show, that, for a generic point p' at the 
MBZ boundary, no symmetry operation relates — p' to 
a vector, equivalent to p'. Hence it is only at point E, 
that the symmetry protects Kramers degeneracy against 
transverse magnetic field. As in the one-dimensional ex- 
ample above, in terms of Eqn. (|5|) this corresponds to 
the simplest case of U = 1 . 

Promotion from rectangular to square symmetry 
brings along invariance under reflections <j\ 2 m cither 
of the two diagonal axes 1 and 2, passing through point 
T in Fig. 0Jb). As a result, eigenstate (Ti#T a U n (7r)|p) at 
momentum 02P (Fig. BJb)) is also degenerate with |p) 
and orthogonal to it, as one can show analogously to the 
examples above. In terms of general condition ([5]). this 
means IA = <y%,%- 

For momentum p at the MBZ boundary in Fig. H^b), 
p and 1T2 P differ by a reciprocal lattice vector and thus 
coincide. Hence, for a square-symmetry lattice in a trans- 
verse field, Kramers degeneracy is protected by symmetry 
at the entire MBZ boundary, as shown in Fig. BJR). 

Degeneracy of special Bloch states in a transverse field 
hinges only on the symmetry of the antiferromagnetic 
state, and thus holds equally in a strongly correlated or a 
weakly coupled material - provided long-range antiferro- 
magnetic order and well-defined electron quasiparticles. 
Under these conditions, quantum fluctuations of the anti- 
ferromagnetic order primarily renormalize the sublattice 
magnetisation, but leave intact the degeneracy of spe- 
cial electron states in a transverse field - certainly in the 
leading order in fluctuations. 

Now, g±(p) can be expanded in a vicinity of the de- 
generacy line g±{p) — 0. With the exception of higher- 
symmetry points, such as point X in Fig. IHR), the lead- 
ing term of the expansion is linear in momentum devia- 
tion Sp from the degeneracy line: 



9±(P) 



Hp • Sp 



(8) 



where 3 p /?i is momentum gradient of gj_(p) at point 
p on the degeneracy line. Inversion symmetry requires, 
that Hp be a pseudo-vector and change sign upon inver- 
sion, and Eqn. ([7]) shows, that S p is of the order of the 
antiferromagnetic coherence length £ ~ 

Zeeman spin-orbit coupling (TTJ) induces a number of 



interesting effects. For instance, substantial momentum 
dependence of g±(p) means, that the Electron Spin Reso- 
nance (ESR) frequency of a carrier in a vicinity of the de- 
generacy line varies along the quasiclassical trajectory. In 
a weakly doped antiferromagnetic insulator, this means 
inherent broadening of the ESR line with doping and, 
eventually, loss of the ESR signal. In fact, this may well 
be the reason behind the long-known 'ESR silence' 12 1 
of the cuprates. Suppression of transverse Pauli suscepti- 
bility is another simple consequence of vanishing g± (p) . 

Finally, momentum dependence of g±{p) allows exci- 
tation of spin flip transitions by AC electric rather than 
magnetic field [131 ] - a vivid effect of great promise for 
controlled spin manipulation, so much sought after in 
spin electronics. Its absorption matrix elements are de- 
fined by S p ~ £ of Eqn. and exceed those of ESR at 
least by two orders of magnitude. According to Eqn. {SJ, 
resonance absorption in this phenomenon shows a non- 
trivial dependence on the orientation of the AC electric 
field with respect to the crystal axes, and on the orienta- 
tion of the DC magnetic field with respect to staggered 
magnetisation. 
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